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Abstract 

A recently developed K-x model for near wall turbulent flows is applied to two severe test cas- 
es. The turbulent flows considered include the incompressible flat plate boundary layer with ad- 
verse pressure gradients and incompressible flow past a backward facing step. Calculations are 
performed for this two-equation model using an anisotropic as well as isotropic eddy-viscosity. 
The model predictions are shown to compare quite favorably with experimental data. 


t Permanent Address: Department of Mechanical Engineering, Stevens Institute of Technology, Hoboken, 
NJ 07030 

§ This research was supported by the National Aeronautics and Space Administration under NASA Con- 
tract No. NASI -18605 while the first and the third authors (ST & CGS) were in residence at the Institute 
for Computer Applications in Science and Engineering (ICASE), NASA Langley Research Center, 
Hampton, VA 23665. 




Introduction 


Two-equation turbulence models have become increasingly popular for the calculation of prac- 
tical aerodynamic flows which can have important applications in the design of advanced aircraft. 
The reason for this popularity is clear: two-equation models incorporate substantially more turbu 
lence physics — and require less ad hoc empiricisms — than the older algebraic eddy-viscosity 
models without most of the added difficulties that arise from the computational implementation of 
second-order closure models. Among the existing two-equation models, the K-e model is probably 
the most popular; its successes and limitations have been fairly well documented [1,2]. One major 
difficulty with the K-e model that has yet to be fully resolved involves its application to near wall 
turbulent flows where inaccuracies and numerical stiffness problems can arise (see Patel, Rodi and 
Scheuerer [3] for an interesting discussion of these problems). Recently, Speziale, Abid and 
Anderson [4] developed a new K-x model for near wall turbulent flows wherein modeled transport 
equations are solved for the turbulent kinetic energy K and the turbulent time scale x. While the 
use of the turbulent time scale had been proposed earlier within the context of two equation models 
(cf., Reynolds [5] and Zeierman and Wolfshtein [6]), none of these previous studies rendered a 
model that could be integrated directly to a solid boundary with the no-slip condition applied. Spe- 
ziale et al. [4] demonstrated that their new K-x model yielded excellent results — and was compu- 
tationally robust — when integrated directly to the wall in an incompressible flat plate boundary 
layer with zero pressure gradient. The K-x model alleviates the problem of the lack of natural 
boundary conditions for the dissipation rate in the K-e model since the turbulent time scale x van- 
ishes identically at a solid boundary. 

The purpose of the present paper is to provide a more comprehensive testing and evaluation of 
the K-x model. Two test cases are chosen to evaluate the model: (i) the incompressible flat plate 
boundary layer with adverse pressure gradients, and (ii) incompressible flow past a backward fac- 
ing step. Unlike in the previous study [4], the K-x model will be solved with an anisotropic as well 
as isotropic eddy viscosity. It has been demonstrated during the past decade that, in some cases, 
the predictions of two-equation models can be enhanced considerably by the use of an anisotropic 
eddy viscosity [7,8]. The results obtained from the K-x model will be documented in detail and 
recommendations will be made for future research. 


The K-x Model 


In the standard form of the K-x model, the Reynolds stress tensor x.. = uj uj (where u' t is the 
fluctuating velocity) is of the form 
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where the eddy viscosity is given by 
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and m,- is the mean velocity, K = ^ u- u- is the turbulent kinetic energy, e = v du- /dx du- /dx- 

is the turbulent dissipation rate, x = K / e is the turbulent time scale, and is a dimensionless con- 
stant that assumes a value of 0.09. In (2),/^ is a wall damping function that is needed for asymp- 
totic consistency [4]; it is taken to be of the form 
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where Re ( = — is the turbulence Reynolds number andy + = yu^N is the usual wall coordinate 
given that u T is shear velocity and v is the kinematic viscosity. In an attached boundary layer, the 
shear velocity is usually taken to be J\/p (where x 0 is the wall shear stress); for more complicat- 
ed turbulent flows with separation, the shear velocity is assumed to be proportional to K 1/2 . The 
wall damping function goes to unity sufficiently far from the wall. In the anisotropic version of 
the K-x model, nonlinear corrections are added to the eddy viscosity model as developed by Spe- 
ziale [8]. This leads to the following Reynolds stress model 
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are, respectively, the mean rate of strain tensor and its frame-indifferent Oldroyd derivative while 
C D is a dimensionless constant that assumes a value of 1.68. Of course, the standard eddy viscosity 
model (1) is recovered in the limit as Cp -» 0. 

The Reynolds stress models (1) and (4) are solved in conjunction with modeled transport equa- 
tions for K and x which take the form [4]: 
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is a wall damping function that is needed for asymptotic consistency. The constants C el and C e2 
assume the values of 1.44 and 1.83, respectively, whereas for their preliminary calculations, Spe- 
ziale et al. [4] chose 


®tl “ °x2 “ °K 1-36 ^ 

with the understanding that these constants could be fine tuned to a non-equivalent set of values 
when a future optimization over a variety of turbulent flows is conducted. When o x j, <r x2 a °d Gr 
are identically equal, this K-x model reduces to an equivalent K-e model which can be simply ob- 
tained by replacing x with K/e in (7) and replacing the x-transport equation (8) with the equivalent 
equation 
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where o e is equal to the common value of a xl , o X 2 , and c K given in (10). 
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Discussion of Results 

The first problem that we will consider is the incompressible flat plate boundary layer with ad- 
verse pressure gradients — the Samuel and Joubert [9] test case. This represents a rather severe 
test since the adverse pressure gradient is strong enough so that the flow is not that far removed 
from separation. In Figure 1(a), the mean velocity profiles predicted by the standard K-x model 
using the eddy-viscosity representation (1) are compared with the experimental data [9] at two sta- 
tions: x = 1.76 m and x = 2.87 m (this corresponds to free-stream Reynolds numbers Re x in the 

range of 2xl0 6 to 4xl0 6 ). It is clear from these figures that the model compares quite favorably 
with the experimental results. In Figure 1(b), the skin friction predicted by the K-x model is com- 
pared with the experimental data. Again, the agreement between the model predictions and the ex- 
perimental data is extremely good. Unlike most other two-equation models, the K-x model does 
not overpredict the skin friction for this adverse pressure gradient test case (see, Rodi and Scheuer- 
er [10] for an interesting discussion of this issue). It should be noted that the use of an anisotropic 
eddy-viscosity does not yield a significant change in the results for this case. 
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The second problem to be considered is incompressible turbulent flow past a backward-facing 
step — the same test case considered by Driver and Seegmiller [11] in their recent experiments. 
For this flow configuration the aspect ratio Hj/H is 8 and the Reynolds number based on the inlet 
free stream velocity and step height is approximately 38,500 (see Figure 2). The calculations were 
conducted on a 200x100 nonuniform mesh using a finite volume method (see Thangam and Hur 
[12]). Here, we use the anisotropic eddy-viscosity model (4) since it has been recently demonstrat- 
ed that such models give rise to improved predictions in recirculating flows [8,12], In Figure 3(a), 
the computed streamlines obtained from the anisotropic K-x model are shown. They indicate re- 
attachment at approximately 6 step heights downstream of separation — a result that is in agree- 
ment with the experimental data [11]. In Figure 3(b), the mean velocity profiles obtained from the 
anisotropic K-x model are compared with the experimental results; again the agreement is extreme- 
ly good. In Figure 4(a), the profiles of the turbulence intensity obtained from the anisotropic K-x 
model are compared with the experimental data. While there are some discrepancies between these 
results in the separation and recovery zone, on balance, the agreement is fairly good. In Figure 
4(b), the profiles of the Reynolds shear stress obtained from the anisotropic K-x model are com- 
pared with the experimental data [11]. Again, while there are some discrepancies in the results 
within the separation and recovery zone, the model predictions are still satisfactory. There is no 
question that this near wall model does a substantially better job than the standard K-e model which 
underpredicts the reattachment point by an amount of the order of 25-30%. 


Conclusion 

The near wall K-x model of Speziale et al. [4] has been tested in two applications involving 
incompressible turbulent flows and the results obtained are quite encouraging. The standard K-x 
model — with an isotropic eddy viscosity — yielded excellent results for the challenging incom- 
pressible adverse pressure gradient case of Samuel -Joubert [9]. For the separated flow considered 
— namely, turbulent flow past a backward facing step — excellent results were obtained from the 
K-x model based on an anisotropic eddy-viscosity. In fact, we believe that these may be the most 
accurate calculations for this backstep problem that have yet been obtained from a two-equation 
model without the ad hoc adjustment of constants. 

Future calculations are planned for supersonic boundary layers with streamline curvature and 
separation (shock-induced or otherwise). For these flows, the K-x model will most probably have 
to be applied with the anisotropic eddy viscosity model and it may become necessary to incorporate 
dilatational effects into the model. After these calculations are completed, and refinements are 
made accordingly, we will have a much better idea of the range of applicability of the K-x model. 
Nonetheless, at this point the prospects appear to be quite promising. 
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Figure 1. Comparison of the results of the K< Model and experiments for the flat 
plate boundary layer with an adverse pressure gradient ( — K*z Model, 
o experiments of Samuel & Joubert [9]) : 

(a) mean velocity, and (b) skin friction coefficient 
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Figure 2. Schematic of the diagram for turbulent flow past a backward facing step 
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Figure 3. Results obtained from the anisotropic K< Model for incompressible 
turbulent flow past a backward facing step with H^H = 8, Re = 38,500. ( — K ~ t 
Model, o experiments of Driver & Seegmiller [11]): 

(a) mean velocity streamlines, and (b) mean velocity profiles. 
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Figure 4. Results obtained from the anisotropic /C-x Model for incompressible 
turbulent flow past a backward facing step with Hj/H = 8, Re = 38,500. ( — K< 

Model, o experiments of Driver & Seegmiller [11]): 

(a) turbulence intensity profiles, and (b) Reynolds stress profiles. 
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